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ABSTRACT 



We perform a search for neutrinos coincident with GRB 080319B — the bright- 
est GRB observed to date — in a ±l,OOOs window. No statistically significant 
coincidences were observed and we thereby obtain an upper limit on the flu- 
ence of neutrino-induced muons from this source. From this we apply reasonable 
assumptions to derive a limit on neutrino fluence from the GRB. 

Subject headings: neutrinos, gamma-rays: bursts 



Introduction 



On 2008 March 19 at 06:12:49 UT (15:12:49 Japan Standard Time) the Swift Burst 
Alert Telescope recorded a gamma-ray bu rst (des ignated GRB 080319B) at equatorial co- 
ordinates (ra = 217°56",dec = +36°18") (||nAs3). This GRB was simultaiieousl y detected 
with the Konus 7-ray detector on board the Wind satellite ( Golenetskii et al. 20081), and sub- 



seque ntly by detect ors operating at lower bandwidths including "P i of the Skv" (jCwiok et al. 



20071 ). TORTORA flMohnari et al.ll2006r ). the Liverpool Telescope (IGuidorzi et al.ll2006f ). the 
Faulkes Telescope North, Gemini-North, the Hob by-Eberly Telescope, and the Very Large 
Telescope. The burst had a duration of tgo > 50 s (IRacusin et al.ll2008bl ) (mak ing it a "long" 
GRB,) with a fluence of 6.23±0.13 x 10"^ ergs cm'^ in the 20 keV-7 MeV band jRacusin et al. 
It is believed that s uch long GRBs may be due to hypernovae from collapsing Wolf- 



2008af ) 

Rayet stars flWooslevI 119931 ) 



This corresponds to an isotropic-equivalent energy release of 
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^ 10^^ ergs at a luminosi ty distance of 1-9 x 1 0^^ cm (IRacusin et al.ll2008al ). It occurred at 
a redshift of z = 0.937 (IRacusin et al. 2008a). This GRB is remarkable as it is the most 



(INaeve 


2008 


). Some GRB models 


rgies ( 


Learned & Mannheim 


2000; 



Meszarosll2006l : lEichler et al.l 119891 ). and so in this article we describe a search for neutrinos 



at Super-Kamiokande coincident with photonic observations of GRB 080319B. 



2. Super-Kamiokande 



Super-Kamiokande is a water Cherenkov detector located within Mt. Ikeno in central 
Japan, under 2, 700 meters water equivalent rock overburden. It has a cylindrical design, 
holds 50 kilotons of water, and is divided into two optically separated sections by a structural 
framework that supports photomultiplier tubes (PMTs). Super-Kamiokande has an inner 
detector (ID) equipped with 11,146 50 cm PMTs aimed inward (at the time considered here,) 
and an outer detector (OD) volume instrumented with 1,885 20 cm PMTs aimed outward 
and equipped with wavelength-shifting plastic plates. The OD functions primarily as a veto 
counter, tagging charged particles that enter or exit the ID. Within the ID we define a 
central 22.5 kiloton fiducial volume within which detector resp onse is expected to be uniform 
( iFukuda. et al. (The Super-Kamiokande Collaboration)! |2003| ). 



Neutrino events (interactions recorded in the detector) with total deposited energy above 
^ 100 MeV are overwhelmingly due to atmospheric neutrinos (decay products from high- 
energy cosmic ray interactions in the upper atmosphere.) The Earth is essentially transparent 
to such neutrinos up to energies on the order of 100 TeV. 

There are several types of events observed at Super-Kamiokande categorized according 
to their energy and/or topology. The low-energy sample is typically used to study solar 
neutrinos and to search for core-collapse supernova neutrinos via elastic u — e scattering for 
all flavors of neutrinos and i^e+P~^ e~^+n for anti-electron- neutrinos. The low-energy event 
selection algorithm — designed to reduce background from cosmic-ray muons, decay electrons, 
penetrating muon-induced s pallat i on events, and radioactivity from the surrounding rock — is 
described in ICravens et al.l (120081 ): iHosaka et al.l (120061 ). Low-energy events associated with 
core-collapse supernovae produced in hypernovae scenarios are not thought to be beamed like 
high-energy neutrinos, making low-energy neutrino detection difficult at these cosmological 
distances. If we assume that rriu = 1 eV and Ei, = 10 MeV, th en the relativi stic delay for 
neutrinos from GRB 080319B would be on the order of 180 s (ILi et al.l l2005l ). which falls 
within the ±l,OOOs search window discussed below. 
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Fully-contained (FC) neutrino events are those where interaction products are observed 
in the ID with no significant correlated activity in the OD, while partially-contained (PC) 
events are those where some interaction products exit the ID. 

Upward-going muon (upmu) events are those where a penetrating particle traveling in 
the upward direction enters and either stops in or passes through the detector, and are 
attributed to muons produced by neutrino interactions in the surrounding rock. We require 
that upmus traverse at least 7 m in the ID if exiting or that they exhibit the equivalent 
energy loss (1.6 GeV) for that range if stopping. Using Monte Carlo data, we measure the 
efficiency of the upmu reduction algorithm by calculating the ratio of events that are both 
true upmus and also tagged as upmus, to events that are true upmus (and may or may 





98% ( 


TI 


iranel 


20 


081). A variety of astrophysical results obtained with upmus can be found 


in 


Fukuda 


(2002 


); 


Desai 


(2008): 


Abe. et al. (The Super-Kamiokande Collaboration) (2006): 


Desai. et al. fThe Super-Kamiokande Collaboration) 


(2004 


)• 



In general terms, FC, PC, and upmu events represent successively higher energy samples 
of neutrino interactions, ranging from ~ 100 MeV and corresponding to > 200 photoelectrons 
(pe) for the lowest energy FC events to above Ei, ~ 1 TeV (< 2 x 10^ pe) for the highest 
energy upmus. Each category of event can be considered a coarse energy bin so long as 
it is understood that the bins overlap significantly. Further details regarding the Super- 
Kamiokande detector design, operation, calibration, a r id da t a redu c tion c an be found in 
Fukuda. et al. (The Super-Kamiokande Collaboration)! (120031 ): lAshid (120051 ). 



The effective area of the detector depends upon zenith angle and varies between 960 — 
1, 300 m^. For each zenith angle, we determine effective area by performing a 2D projection 
of the cylindrical detector onto a plane, subject to the constraint that a particle passing 
through the detector must be able to traverse at least 7 m through the ID. 



3. Search Method 



Following a method developed in a previously published GRB search (IFukudal |2002| ) . 
we employ ±l,OOOs timing window centered on the beginning of photonic observations of 
GRB 080319B and look for coincident neutrino events. This window size allows for any 
reasonable delay, positive or negative, between neutrino emission and photonic emission, 
while still providing a very small likelihood of random coincidences due to atmospheric 
neutrinos. We consider all four categories of atmospheric neutrino events — low-energy, FC, 
PC, and upmu — but the strongest hmits come from the upmu dataset. 
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To eliminate cosmic ray background, we require upmu events to come from below the 
horizon, which means that the upmu sample is only sensitive to declinations less than +54°. 
Fortunately, GRB 080319B occurred 17° below the horizon, with a local azimuthal angle of 
= 32°. 

The upmu reduction algorithm is divided into several stages design ed to filter out 
downward-going cosmic-ray muons, which sometimes masquerade as upmus (lAshiell2005l ). as 
well as ultra-high-energy (UHE) events, (with > 1.75 x 10^ pe in the ID,) that are analyzed 
separately for astronomy studies. UHE background events occur at a rate of about 30 day~^. 
Upward-going UHE events, ( which can be att ributed to neutrinos,) are very rare and occur 
at a rate of 2.7 x 10~''day~^ (ISwanson 1 120061 ). 



Currently favored GRB models predict a hard d^i,/ dE (x E spectrum (ILearned &: Mannheim 



2000l : iMeszarod |2006| ) . Combining this with the roughly linear energy dependence of the 
neutrino-nucleon cross section, plus the fact that the long range of high-energy upmus makes 
the effective detector volume increase with energy, we expect a peaked spectrum for observed 
muons, as shown in Figure [H The position of the peak depends upon how hard the GRB 
neutrino spectrum turns out to be, but the most likely range for spectral indices (7 = 2 ~ 3) 
puts the peak in the energy region spanned by Super-Kamiokande's upmu events. 

The detector's angular resolution — defined such that 68% of events have an angular 
separation between true and fit muon direction that is smaller than the resolution — is about 
1° at this zenith angle. At energies above 1 Te V the typical angular separa tion between a 
neutrino and its daughter muon is less than 1° ( iLearned fc Mannheimll2000l ). The effective 
area at this zenith angle is ~ 1, 270m^. 

The rates of post-reduction atmospheric and solar background events are summarized 
in Table [H Since the background rates are so low, a coincidence in time provides a strong 
signal of a correlation between the observed neutrino and the photonic detection. (The 
probability of an accidental upmu coincidence from an atmospheric neutrino is about 4%.) 
It is also possible to employ a cut on the reconstructed upmu direction in relation to the 
GRB direction, but this was not necessary as there were no upmu coincidences in time. 



4. Results 

No statistically significant coincidences were observed in any of the four categories (see 
Table [1]). Therefore we set an upper limit on neutrino-induced muon fiuence, the effective 
neutrino flux integrated over the duration of the burst, which is the appropriate physical 
quantity to describe transient sources. We find that the upper limit at 90% on the fiuence of 
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neutrino-induced muons from GRB 080319B is <l>f'^'' = 1.96 ± 0.04 x 10"^ cm-^. (The ^ 2% 
uncertainty is due primarily to uncertainty in our track length reconstruction algorithm and 
in our live time calculation.) 



5. Neutrino Fluence Limit 

Extrapolating from upmu fluence limits to neutrino fluence limits requires making an 
assumption about the nature of the source spectrum. We estimate the muon neutrino flux 
at the detector (after neutrinos have oscillated from their source distribution.) Since a wide 
variety of models predict a power law spectrum with a spectral index of about 7 = 2, we 
assume d^^/dEi, oc E~'^ and consider the cases of 7 = 2 and 7 = 8. 

Upmu fluence can be related to neutrino fluence as in Equation [H 



dE„ 



JTJrnm 



P{E,,Ef^)S{z,E, 



d^u 
dE, 



Here $^(> E™'^ 



= 1.6 GeV. The quantity P{E,, E"^"" 



is the fluence of upmus with energies above the minimum upmu energy of 

is the probability that a neutrino with energy 
El, creates a muon with energy greater than -E™" and S{z, E,) is the Earth's shadow fa ctor. 
We use cross sections from the GRV 94 parton distribution f unction (|Gluck et al.lll995l ) and 
the muon range is determined using iLipari &: StanevI (Il99ll ): iRend (120051 ). 



) 



PiEu, E^ 
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Na I dE^ 
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Table 1. Background rates and events observed in window. 



type 


background (day ^) 


...in itljOOOs window 


observed in window 


low-energy 


390 


9 


11 ±3 


FC 


8 


0.2 





PC 


0.6 


0.01 





upmu 


1.6 


0.04 





UHE 


2.7 X 10"* 


6 X 10"'' 
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d<I>u/dE (arbitrary units) vs. log(E) 




2 4 6 8 



Log(E[GeV]) 

Fig. 1. — The upmu spectra for a source with spectral indices 7 = 2 (dashed) and 7 = 8 
(sohd). Analysis cuts impose an energy threshold E^™ > 1.6 GeV. 
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Here dace /dE^{E^^ E^) is the charged current component of the neutrino-nucleon cross 
section; (neutral current interactions do not produce penetrating muons.) The quantity 
R{Ef^, -E™™), meanwhile, is the average range in rock for a muon with an initial energy of 
and a final energy greater than i^;™™. The quantity is the water equivalent Avogadro's 
number (scaled by the density of water.) 



S{z, By) = e" 



(3) 



Here lcoi{z) i s the angle-dependen t column depth along the neutrino path. Column depth is 



calculated in Gandhi et al. 



(119961 ) using the "Preliminary Earth Model." The results of this 



calculation are summarized in Table [2J 

In order to derive fluence limits for low-energy and FC/PC events we use Equation |U 
which differs from Equation [1] because low-energy and FC/PC interactions must occur in 
the water fiducial volume, while upmu interactions occur in the rock outside the detector. 



FC/PC 



90 



(4) 



Here A^^go is the 90% CL limit on the total number of FC/PC neutrino interactions in the 
window and is the number of interaction targets of type i where i can designate neutrons, 
protons, or electrons. The quantity a'^{E^) is the cross section for neutrinos to interact with 
particles of type i and e{E,^) is the detec tor efficiency as a function of neutrino energy. The 
low-energy limits use cross sections from lStrumia fc Vissanil (120031 ). The FC/PC limits and 
low-energy limits are presented in Table [31 



6. Systematic Uncertainty 

Uncertainty in neutrino-nucleon cross section creates systematic uncertainty in the cal- 
culated neutrino fluence. In Equation [H we utilize a table of neutrino- nucleon cross section 



which has an associated uncertainty of ~ 10% at the upmu energy scale (iGandhi et al.lll996l ). 



Table 2. Limits at 90% CL on the fluence of muon neutrinos and antineutrinos from 
GRB 080319B given different spectral indices and derived with the upmu sample. 



spectral index (cm ^) (cm ^) 



7 = 2 
7 = 3 



16 ± 1.7 22 ±2.4 

6.6 ±0.7x10^ 10±l.lxl03 
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To assess how this affects the neutrino ffuence uncertainty, we carry out the calculation us- 
ing $_|_ = $(o" —>■ 110% cr) as well as $_ = $(cr —>■ 90% cr) and estimate the uncertainty as 
half the difference. For upmus, an additional uncertainty arises from cross section uncer- 
tainty in the Earth shadow. We add these two uncertainties in quadrature and find the total 
uncertainty in $y to be 11%. 



7. Comparison with Other Results 



In iFukudal (120021 ) a set of 1,454 GRBs from the BATSE catalog was used to test for 
coincident neutrinos at Super-Kamiokande. Assuming a spectral index of 7 = 2, a 90% 
limit was obtained for the average GRB fiuence of $^0% < o.038 cm ^ for (0.050 cm ^ 
for z/^.) Since this limit is derived from the summed fiuence of 1,454 GRBs it is, of course, 
significantly lower than the one quoted here and obtained from the consideration of a single 
GRB. (The lim it quoted here, however, is better than the "average" single-burst limit from 
Fukudal J2002h . $^ < 55 cm-^.) 



All the same, the importance of the present result is magnified by the fact that GRB 080319B 
was so exceptionally bright. Further, studies of burst-to-burst fluctuations suggest that vari- 
ations in the Lorentz factor characteri zing GRB shockwayes ca n change the neutrino fiuence 
by roughly four orders of magnitude (IHalzen fc Hooperlll999l ). L ikewise, modeling o f indi- 
vidual GRBs reveals a similarly wide range of predicted fiuences ( iGuetta et al.ll2004l ). 



The best limits on average GRB fiuence currently come from the AMANDA experiment, 
which reports an upper limit of 1.4 x 1 0~^ cm~^ for neutrino en ergies between 250 — 10^ GeV 



and assuming a spectral index of 7 = 2 (lAchterberg et al.ll2008l ) . AMANDA has also searched 
for neutrinos from a single GRB (GRB 030329) and se t differential fiuence limits ( for several 
models) in the range of d^/dlogE = 1 — 6 GeVcm"^ (jStamatikos &: Clarkd 120051 ) . 



8. Conclusions 

We have looked for neutrinos coincident in time with GRB 0803 19B — the brightest GRB 
observed to date. We found no significant coincidences, and so we set an upper limit on the 
fiuence of neutrinos from this object. 
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Table 3. Limits at 90% CL on the fluence of neutrinos from GRB 0803 19B given a 
spectral index of 7 = 2 and derived with the FC/PC data sample (top), with the 
low-energy sample (middle) , and with UHE data (bottom) . 



species (cm ^) 

from FC/PC only 

i^^ 3.5 X 10"' 

7.6 X 10" 

Ue 4.3 X 10-* 

Pe 7.6 X 10* 

from low-energy only 

5.1 X lO'' 

Ve 8.0 X 10'^ 

Pe 3.0 X 10* 

from UHE only 
u^, 1.2 X 103 

u^, 1.6 X 10^ 



